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Abstract 

The behaviour of the transverse momentum fluctuations with the cen- 
trality of the collision shown by the Relativistic Heavy Ion Collider data is 
naturally explained by the clustering of color sources. In this framework, 
elementary color sources -strings- overlap forming clusters, so the number 
of effective sources is modified. These clusters decay into particles with 
mean transverse momentum that depends on the number of elementary 
sources that conform each cluster, and the area occupied by the cluster. 
The transverse momentum fluctuations in this approach correspond to the 
fluctuations of the transverse momentum of these clusters, and they behave 
essentially as the number of effective sources. 
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Event-by-event fluctuations of the mean transverse momentum are considered to be 
one of the most important tools to identify a phase transition in the evolution of the 
system created in relativistic heavy ion collisions, since a second order phase transition 
may lead to a divergence of the specific heat which could be observed as fluctuations in 
mean transverse momentum [1]- [3]. These fluctuations have been extensively studied 
both theoretically [4]- [10] and experimentally [11]- [17]. Recently, the PHENIX collab- 
oration [15]- [16] of the Relativistic Heavy Ion Collider (RHIC) has reported a peculiar 
behaviour of the variable Fp^ that measures the transverse momentum fluctuations 
as a function of the number of participants in Au-Au collisions. F^^ quantifies the 
deviation of the observed fluctuations from statistically independent particle emission, 

J-, ^data ^random / -, \ 

^PT = 

'■^random 

where 



'<pI>-<Pt>^ 

u = , (2) 
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and < pt > is the mean transverse momentum averaged over all particles and all 
events. The data show that Fp^ increases with the number of participants, reaching a 
maximum around Npart — 150 -j- 200 and decreasing at higher centrality. In this paper 
we show that this behaviour is naturally explained by the clustering of elementary color 
sources that may take place in heavy ion collisions. 

Multiparticle production is currently described in terms of color strings stretched 
between the partons of the projectile and the target. These strings decay into new 
ones by sea q — q production, and subsequently hadronize to produce the observed 
hadrons. In our approach, the strings are equivalent to effective color sources with 
a flxed tranverse size vrrQ, with ro ~ 0.2 fm, fllled with the color held created by the 
colliding partons. With increasing energy and/or atomic number of the colliding nuclei. 



the number of exchanged strings grows, so they start to interact forming clusters. In 
the transverse space that means that the transverse areas of the strings overlap, as it 
happens for disks in the two dimensional percolation theory. Moreover, at a certain 
critical density of disks, rj ~ 1.12 — 1.18, a macroscopical cluster appears which marks 
the percolation phase transition [18]- [19], which is a second order, non thermal, phase 
transition. 

In the case of a nuclear collision, the density of disks -elementary strings- corre- 
sponds to 

NsSi 

Vc = (3) 

where Ng is the total number of strings created in the collision, each one of an area 
= Trro, and Sa corresponds to the nuclear overlap area, Sa — t^R\ for central 
collisions. 

The percolation theory governs the geometrical pattern of the string clustering. Its 
observable implications, however, required the introduction of some dynamics in order 
to describe the behaviour of the cluster formed by several overlapping strings [20]- [21]. 
We assume that a cluster of n strings that occupies an area Sn behaves as a single 
color source with a higher color field, generated by a higher color charge Qn- This 
charge corresponds to the vectorial sum of the color charges of each individual string 
Qi. The resulting color field covers the area Sn of the cluster. As = (Si Qi)^, 
and the individual string colors may be oriented in an arbitrary manner respective to 
one another, the average QijQij is zero, so — nQ\. Qn depends also on the area 
Si of each individual string that comes into the cluster, as well as on the total area 
of the cluster Sn-, Qn — \f^§^Qi *■ We take constant and equal to a disk of radius 



"Qn would be equal to ^/nQ\ if the strings overlap completely. Since the strings may overlap 



To ~ 0.2 fm. Sn corresponds to the total area occupied by n disks [20]. One could do 
reasonable alternative assumptions about the interaction among the strings, but they 
have incompatibilities with correlation data [22]- [23]. 

Notice that if the strings arc just touching each other, Sn = nSi and Qn = nQi, so 
the strings behave independently. On the contrary, if they fully overlap, Sn — Si and 
Qn = y/nQi- Knowing the color charge Qn, one can compute the multiplicity /i„ and 
the mean transverse momentum < >„ of the particles produced by a cluster of n 
strings. According to the Schwinger mechanism for the fragmentation of the cluster, 
one finds 

InSn , /n5'iNi/4 

for the multiplicity /i„ and the average transverse momentum < of the particles 

produced by a cluster formed by n strings, where /xi and < pr >i correspond to the 
mean multiplicity and the mean transverse momentum of the particles produced by 
one individual string. These equations constitute the main tool of our evaluations. 

The behaviour of the transverse momentum fluctuations can be understood as fol- 
lows: At low density, most of the particles are produced by individual strings with the 
same < pr >i, so the fluctuations are small. Similarly, at large density above the per- 
colation critical point, there is essentially only one cluster formed by most of the strings 
created in the collision and therefore fluctuations are not expected either. Instead, the 
fluctuations are expected to be maximal below the percolation critical density, where 
the number of clusters is larger. Moreover, there are clusters formed by very different 
numbers of strings, with different size, and therefore with different < pt >n- 

only partially we introduce a dependence on the area of the cluster. See first Ref. of [20] for 
more details. 
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In order to develop quantitatively this idea, we introduce the function [4] (f) defined 

by 



'< Z2 > 



V < A* > 

Fp^ is related to (f) [14], approximately 



- \/< > . (5) 



For each particle we define Zi — Pn— < Pt >, where pxi is the transverse momentum 
of the particle i and < pt > is the mean transverse momentum of all particles aver- 
aged over all events. \J < > is the second moment of the single particle inclusive z 
distribution, and it is averaged over all events. Z is defined for each event, 

Z^^Y.h (7) 

where Ni is the number of particles produced in an event i. 

In this way, introducing our formulae for the multiplicity and the mean pt we get: 

< » EjlT"'" Ej l^Uj < Pt >nj .gs 

The sum over j goes over all individual clusters j, each one formed by Uj strings and 
occupying an area The quantities rij and are obtained for each event, using 
a Monte Carlo code [24]- [25], based on the quark gluon string model. Each string 
is generated at an identified impact parameter in the transverse space. Knowing the 
tranverse area of each string, we identified all the clusters formed in each event, the 
number of strings rij that conforms each cluster j, and the area occupied by each 
cluster Snj- Note that for two different clusters, j and k, formed by the same number 
of strings rij — Uk, the areas Snj and Sn^. can vary. Because of this we do the sum over 
all individual clusters. 



5 



For the quantities < pr >nj -mean px of the particles produced by a cluster j of 
rij strings and area Sn^- and -mean multiplicity of a cluster j formed by rij strings 
and of area Snj- we apply the analytical expressions given by eqs. (4). Finally we do 
the average over all events. 

By introducing eqs. (4) for < px >nj and we get: 



< Pt >= 



n^Sn, \l/2 
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where the mean value in the r.h.s. corresponds to an average over all events. 
For the quantities < z'^ > and < > we obtain: 

Efir- E, (^) "\.[m "'<Pr>,-<PT> 
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and 
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Hi < Pt >i 
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Finally we arrive to: 
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Pt ~ 




h + /4/I - 2/2/5 
/i - /I 



- 1 . 



(12) 



where /i, /2, /a, /4 and /s are defined in the expressions (9)- (11). 

In order to compute eq. (12), several ingredients are necessaries. On one hand we 
need a Monte Carlo code [24]- [25] for the cluster formation, in order to compute the 
number of strings that come into each cluster and the area of the cluster. On the other 
hand, we do not use a Monte Carlo code for the decay of the cluster, since we apply 
analytical expressions (eqs. (4)) for the transverse momentum and the multiplicities of 
the clusters. 

We also need the value of ii\ -multiplicity produced by one individual string-. It 
was previously fixed from a comparison of the model to SPS and RHIC data [20]- 
[21] on multiplicities. In the first Ref. of [20], the total multiphcity per unit rapidity 
produced by one string has been taken as /xotot — 1- If we assume that 2/3 of the created 
particles are charged, that would lead to a charged particle multiplicity per unit rapidity 
for each individual string of /xq c/i = 0.65. In order to compare with experimental data 
we define //i = /J'OchV, where y is the rapidity interval of the produced particles. We 
don't introduce any dependence of Hq with the energy or the centrality of the collision. 
Notice that in (12) the value of < >i cancels. 

We have neglected the subsequent rescattering of hadrons and resonances that takes 
place after the decay of the clusters. It gives rise to correlations which would be similar 
in the clustering approach and in an independent string picture, unless additional 
dynamics were taken into account. Therefore its contribution to Fp^ cancels in our 
approach. 

The comparison of our results for the dependence of Fp^, on the number of partici- 
pants Np with the PHENIX data [16] is shown in Fig. 1. The calculation is done for 
charged particles in the rapidity range I77I < 0.35, /ii — 0.7 /loch- An acceptable overall 
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agreement is obtained. 

In order to compute our value for Fp^,, we take into account all possible transverse 
momenta, whereas in the experiment there is a limited acceptance, 0.2 GeV/c < pr < 
pTnax^ PHENIX [15]- [16] has studied the variation of Fp^ with the maximal value of the 
acceptance for px, p^:"^. The maximum of Fp^ is reached for the largest acceptance, 
p-max _ ^ QgV/c [15]. So we Can expect that our value for Fp^ is going to be higher 
than the experimental one, specially for a moderate number of participants, Np, since 



number of participants for pj^*" > 2 GeV/c. This means that, for momenta higher than 
2 GeV/c, the high px contribution would be due to collisions with a moderate number 
of participants. These considerations may explain the difference between our results 
and PHENIX data -with a limited acceptance of 0.2 GeV/c < pt < 2.0 GeV/c- at low 



In Fig. 2 our results for (pp^ of charged particles in Pb-Pb central collisions at 158 
AGeV are compared with the experimental data of NA49 Collaboration [27]. In this 
case the data correspond to the forward rapidity range 4.0 < ?/ < 5.5. For this reason 
we use //i = l.h ^Qch-i which in principal implicates larger correlations. However we see 
that this effect is compensated, since we have a lower value for the mean number of 
strings at fixed Np, due to: 

a) lower energy at SPS than at RHIC so less strings are produced, 

b) the mean number of strings in this rapidity region is proportional to the number 
of participants Np, while in the central region it is proportional to N^^^ due to the 
contribution oi q — q strings from the sea. 

For the computation of 0^^, we use < pt >i— 0.3 GeV/c. 

The CERES Collaboration has also measured 4>p^ [27] at four different centralities: 
- 5%, 5 - 10%, 10 - 15% and 15 - 20% for Pb-Au coUisions at 40, 80 and 158 
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the truncated average pt [26], < > 



/„min PTdN/dpT 



/„min dN/dpT 




— p™*", decreases with the 



AGeV/c in the central pseudorapidity range 2.2 < 77 < 2.7 and restricted to tracks 
with transverse momenta 0.1 < < 1.5 GeV/c. Due to the narrower rapidity and 
transverse momentum range one would expect a lower (pp^ value. However this is 
compensated with the increase of the number of strings at central rapidities. The 
data at 158 AGeV, after short range removal, are 3.3, 3.6, 4.4 and 4.1 for the above 
mentioned centralities, with errors of the order of 1.5 -for smaller energies the data are 
lower as expected-. These values and their dependence with centrality are compatible 
with our results of Fig. 2. 

In order to have a better understanding of the behaviour of Fp^, and (pp^. on the 
number of participants, wc plot in Fig. 3 and Fig. 4 the mean number of clusters 
M and the dispersion on the number of clusters multiplied by the number of clusters 
(7m * M at RHIC and SPS energies. The ratio (jIj/M would be one in the case of a 
Poisson distribution. The pr fluctuations are due in our approach to the different mean 
transverse momenta of the clusters. These momenta depend on the number of strings 
that comes into the cluster and the area occupied by the cluster through our eq. (4), 
therefore M and gm should be the key quantities. However, a\j/M ranges between 
1/2 and 2 in the whole Np range, what indicates a lower variation than the one for M, 
as can be seen from Fig. 3 and 4 where M and gm * M are plotted. The only effect 
of (7m is to shift the maximum of M. Because of this we expect the dependence of 
Fp^ and 0^^, on Np to be more similar to the M behaviour, as it is actually. In other 
words, a decrease in the number of effective sources leads to a decrease of the tranverse 
momentum fluctuations. 

Similar conclusions have been reached in Ref. [29], where a formula for Fp^ as a 
function of the cluster dispersion over the mean number of strings per cluster has been 
obtained. 

Notice that we only need to know the number of strings formed for each centrality 
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and their location in the impact parameter space in order to form clusters. This 
information, together with eq. (4), is enough for us to calculate F^^, and 0^^,. The 
same variables have been able to describe the behaviour of the strength of two [22] 
and three [23] body Bose-Einstein correlations with centrality and the dependence of 
the multiplicities and transverse momentum distributions [21], [28] on the centrality. 
All that points out that the percolation approach may be appropriate to describe the 
relativistic heavy ion collisions. 

We thank J. Dias de Deus for useful discussions. This work has been done under 
Contract No FPA2002-01161 from CICYT of Spain and FEDER from EU. 
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FIGURES 




50 100 150 200 250 300 350 400 



Number of participants 

FIG. 1. Fpj,{%) versus the number of participants. Experimental data from PHENIX at 
= 200 GeV are compared with our results (solid line). 
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FIG. 2. (pp^ versus the number of participants. Experimental data from NA49 Collabo- 
ration at SPS energies are compared with our results (solid line). 
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FIG. 3. Mean number of clusters M versus the number of participants for Pb-Pb collisions 
at SPS energies (dotted line) and Au-Au collisions at RHIC energies (solid line). 
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FIG. 4. Dispersion on the number of clusters multiplied by the number of clusters aM*M 
versus the number of participants for Pb-Pb collisions at SPS energies (dotted line) and Au-Au 
collisions at RHIC energies (solid line). 
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